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BACKGROUND OF THE INVENTION 

FIFJ D OF TI IR INVENTION Field of the Invention 

The present invention relates to a method of evaluating a capillary pressure curve of 
the rocks of an underground reservoir from measurements on rock debris taken 
5 therefrom. 

RAKfiROT INF) OF THE TNVENTlON Description of the Prior Art 
Laboratory measurements on cores or cuttings 

Measurement of petrophysical parameters such as the permeability, the porosity and 
the capillary properties on rock fragments obtained while drilling a well through an 
10 underground formation is an interesting opportunity for operator companies to rapidly 
obtain a first petrophysical characterization of producing zones traversed by the well. 

French Patent FR-2,809,821 i filed by the applicant assignee, describes a system of 
evaluating physical parameters such as the absolute permeability of porous rocks of an 
underground reservoir zone, from cuttings returned to the surface in the drilling mud. In 

15 an enclosure where the cuttings are dipped in a viscous fluid, some of this fluid is 
injected at a pressure that increases with time until a predetermined pressure threshold is 
reached, so as to compress the gas trapped in the pores of the rock. This injection stage 
is followed by a relaxation stage where injection is stopped. The pressure evolution 
during the injection process being-is modelled from initial values selected for the 

20 physical parameters of the cuttingsr-a- . A computer adjusts them-the values iteratively 
so as to obtain the best possible agreement between the modelled pressure curve and the 
pressure curve really measured. 
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French P atent- patent application FR-02/0023 a filed by the applicant assignee. 
describes another method of evaluating physical parameters such as the absolute 
permeability and the porosity of the rocks of an underground reservoir zone, also from 
cuttings. An enclosure containing the rock fragments and filled with a viscous fluid is 

5 communicated communicates with a vessel containing the same fluid at a predetermined 
pressure so as to compress the gas trapped in the pores of the rock. The time of 
application of this pressure, according to whether it is short or long, allows to measure 
measuring either of the pressure variation in the enclosure or the variation of the volume 
actually absorbed by the rock fragments. The pressure or volume evolution in the 

10 enclosure is then modelled from initial values selected for the physical parameters of the 
fragments, and the values of the physical parameters of the rock fragments are 
iteratively adjusted so that the modelled evolution best adjusts to the measured 
evolution of the physical parameter in the enclosure. 

In the sphere- field of petrophysical charact e rization study , the capillary pressure is 
15 also a very important datum for operators because it conditions : 

- the initial distribution of the fluids in the reservoir from the aquifer zone (referred to 
as WOC, for water-oil contact, by the man those skilled in the art) to the upper part of 
the reservoir (transition zone). According to the capillary pressure curve associated with 
a reservoir rock and the nature of the fluids in place, this transition zone can extend over 
20 some — many meters or s ome ten meters , which has an important effect on the 
determination of the accumulations in place, 

the input pressure of a rock, which is particularly important for cap rocks. For 
example, for a gas storage tank, the input pressure of the cap rocks conditions directly 
the allowable overpressure in the storage levels without having leaks. 
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With the current techniques, the capillary pressure curve is obtained by means of 
laboratory measurements on reservoir cores. These methods are expensive because of 
the coring operations as well as the measurements on the cores, and the results are often 
available only several months after drilling. 

5 Approaches to rapidly obtain the capillary pressure curve 

However, there are alternative methods described in the literature for evaluating the 
capillary pressure curve rapidly, either during drilling or slightly later. 

The most commonly used approach consists in using uses the mercury porosimetry 
technique for measuring the air/mercury capillary pressure curve Pc directly from 
10 cuttings. However, the obtained curve obtained is significantly different from the 
reference curve obtained on— from cores with high wetting fluid saturations. 
Besides Further , this approach is based on the use of mercury, which is extremely 
polluting and progressively forbidden by the law in many countries, which poses a 
major problem for applying this technique in the near future. 

15 Another known method uses the Nuclear Magnetic Resonance (NMR) technique to 

rapidly estimate the capillary p ressure curve from strati graphic d ata measured i n t he 
well shortly after drilling. It is notably described in the following publications : 

- Bowers, M., A. et al. : « Pr e diction of permeability from capillary pr e ssure curves 
derived with NMR » "Prediction Of Permeability From Capillary Pressure Curves 

20 Derived With NMR" , 17 September 1998, 

- Marshall, D., et al. : « Method for correlating NMR r c laxom c try and m e rcury 
inj e ction data » "Method For Correlating NMR Relaxometry And Mercury Injection 
Data" , SCA No. Society of Core Analysts International Symposium 1995, 
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- Volokitin, Y., W. J. et al. : « A p ra ct i c al approach to obtain 1 st drainage capillary 
pressure ewves from NMR cor e and log data » "A Practical Approach To Obtain 1 st 
Drainage Capillary Pressure Curves From NMR Core and Log Data" , SCA No. Society 
of Core Analysts International Symposium 1999. 

5 The NMR relaxation signal is first converted in terms of pore size distribution, then 

in terms of threshold size distribution, which allows to calculate calculation of a 
pseudo-capillary pressure curve. This approach has been tested on several samples of 
known curve Pc. The results show that a good agreement with the reference curves can 
only be obtained by means of a rigorous calibration stage to be carried out case by case 

10 according to the nature of the rocks studied. This calibration stage is necessary owing to 
the uncertainty on : 

- the NMR signal-pore size distribution conversion which depends on the value of the 
surface relaxivity which is variable according to the rocks, and 

- the pore size distribution-threshold size distribution conversion which depends on 
15 the nature of the rock and on the diagenesis process. 

This approach is therefore not recommended in a predictive exploration context. In 
any case, it would not be applicable to cuttings. 

Image analysis has also been the subject of work intended to obtain a curve Pc. The 
porous medium is first prepared in a_form of a thin section photographed by a_scanning 
20 electron microscopy or SEM. The image obtained is then analys e d analyzed so as to 
determine parameters representative of the proportion and of the shape of the voids in 
relation to the rock. In particular, it is possible to determine a threshold size distribution 
to reconstruct a pseudo-capillary pressure curve Pc. The main limitation of this method 



is the two-dimensional (2D) nature of the thin section, whereas the capillary pressure is 
by definition a three-dimensional (3D) property. Besides, this technique requires quite 
heavy conditioning, which is not really compatible with a result obtained slightly later. 
Image analysis could be applicable en-to cuttings but it would require careful calibration 
5 to acquire a good predictability. 

Finally, it can be noted that the centrifuging technique is sometimes applied in the 
field to cuttings , but in ord e r in order to extract the largest possible amount of drilling 
fluid from the rock to minimize pollutant discharges to the environment and to limit the 
cost by recycling the drilling fluid recovered. A s far as wo know, no The inventors do 
10 not know of centrifugation en— of cuttings has -having b een considered in order to 
determine capillary properties. 

SUMMARY OF THE INVENTION 

The method according to the invention allows to determin e determination of the 
capillary pressure curve of rocks of an underground reservoir from measurements on 
15 rock debris or fragments (such as cuttings) taken therefrom, over the total saturation 
range of these rocks, within a short period and at a low cost, from these measurements. 
It comprises : 

measuring the permeability k of the rock debris^; 

measuring the capillary pressure curve Pc as a function of the saturation of the rock 
20 debris initially saturated with a fluid by subjecting them to centrifugation 7 i and 

- parametrizing a capillary pressure curve Pc satisfying empirical relations depending 
on adjustable parameters, that is constrained to adjust to an asymptotic part of the 



capillary curve measured by centrifugation, and to the value of permeability k measured 
on the cuttings, so as to obtain the-awhole of the capillary pressure curve. 

Parametrizing the curve is advantageously carried out by selecting by default a set 
0 f sakt- the parameters allowing calibration on the asymptotic part of the capillary 
5 pressure Pc with low saturations, and by modifying the parameters step by step so that 
the estimation of the permeability given by one of the empirical relations used is best 
adjusted with the measurements of permeability k carried out on rock debris and-with 
th4s- the asymptotic part. 

Permeability k of the cuttings is measured for example from measurements of the 
10 pressure variations in a vessel filled with a fluid containing the cuttings after it has been 
communicated coupled for a predetermined period of time with a tank containing the 
same fluid under pressure, and from the volume actually absorbed by the cuttings, and 
from modelling the evolution of the pressure or of the volume in the vessel, from initial 
values selected for the physical parameters of the cuttings, which are iteratively adjusted 
15 so that the modelled pressure evolution best adjusts with the measured evolution of the 
physical parameters of the cuttings. 

The method is notably advantageous in that it provides the capillary pressure of the 
rocks on the basis of simple cuttings that are more readily available and less expensive 
to obtain. The results are also obtained much more rapidly than with cores. 

20 BRtEF nRSPRIPTION OF THE F1GURES BR1EF DESCRIPTION OF THE 

DRAWINGS 
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Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of an embodiment given by way o f non 
limitative example, with reference to the accompanying drawings wherein : 

- Figure Fig. 1 shows the evolution of the NMR signal during various centrifugation 
5 stages carried out on cuttings^; 

. Figure 2a Fig. 2A shows the curves Pc obtained on a single rock (B7) from cuttings 
and from a core?; 

- Figure 2b Fig. 2B shows the comparison between the reference curve Pc measured 
on a core and the reconstructed curve Pc over the total saturation range from the 

10 centrifugation measurements and the permeability measurements 

- Figures 3a and 3b Figs. 3A and 3B show results comparable to those shown in 
Figs.2a and 2b Figs. 2A and 2B respectively for another rock GDVU; 

. Figures 1a and 1 b Figs. 4A and 4B show results comparable to those shown in 
Figs.2a and 2b Figs. 2A and 2B respectively for another rock Rotl ? ; 

15 - Figures 5a and 5b Figs. 5A and 5B show results comparable to those shown in 
Figs.2a and 2b Figs. 2A and 2B respectively for another rock St Max ? ; 

- Figure Fig. 6 shows the comparison between the reference permeabilities and the 
permeabilities measured on cuttings within the scope of the method described in the 
aforementioned French patent application FR-02/0023 T ; and 

20 - Figure Fig. 7 diagrammatically shows the flowchart for implementing the method: 
OT7TATT ,FD DHSCRIFTION DETAILED DESCRIPTION OF THE INVENTION 
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The method intended for fast evaluation of a capillary pressure curve Pc from rock 
fragments or cuttings according to the invention is illustrated by Figure Fig. 7. The 
method is based on two e xperimental measuring s tages, followed by a parametrizing 
stage by reference with known curves, to reconstruct curve Pc over the total saturation 

5 range. During one of the measuring stages, capillary pressure data Pc are acquired by 
centrifuging the initially water-saturated rock fragments. The other experimental stage 
allows to calculat e c alculation of the permeability value k of the rock from the method 
described in the aforementioned French patent application FR-02/0023. Reconstruction 
of curve Pc over the total saturation range is carried out by using a parametrized form. 

10 The parameters of the curve are determined in such a way that the curve matches the 
capillary pressure data Pc obtained experimentally and the estimated permeability value 
from the curve Pc obtained according to a known method such as Thomeer's or 
Swanson-Kamath's method, which are remind e d described b elow, with the value 
measured on the rock fragments. Several application cases are presented, with show the 

1 5 very good agreement obtained with reference c urves without any particular previous 
calibration procedure. 

I) Measurement of Pc by centrifugation from the cuttings 

To implement the method, it is possible to use standard centrifugation m eans or 
more sophisticated means with automatic monitoring of the volumes of fluid produced, 
20 such as those described for example in patents EP-603,040 (TJS~ Patent 5,463,894), FR- 
French patents 2,763,690, FR-2,772,477 (US- patent 6,185,985) or FR 2,798,734 filed 
by the applican t assignee . 

The cuttings which have come up to the surface during the drilling operation are 
first cleaned with solvents in a Soxhlet type device, then dried and saturated with 30 g/1 
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brine. The cuttings are then drained in a damp cloth so as to remove the water trapped 
between the various cuttings d, then placed in a cell or cup fastened to the end of a 
rotating arm. The water expelled by centrifugation from the cuttings flows through a 
grate and it is collected at the base of the cup. The experimental data are acquired in the 
5 same way as in the context of centrifugation on a core. For a centrifugation stage (given 
rotating speed), the evolution of the water production is measured until no significant 
variation can be observed any longer, then the rotating speed is increased to start a new 
stage. 

As can be seen in Figure Fig. 1, a progressive decrease of the NMR signal and a 
10 shift to the short relaxation times T2 are observed, which express a desaturation of the 
porous medium with the increase of the rotating speed. This type of measurement shows 
that there definitely is a capillary contact between the cuttings, which allows to measure 
Pc by centrifugation from cuttings. 

The volume of water produced during the experiment is converted to saturation data 
15 from the volume of water initially contained in the cuttings. The latter is determined by 
weighing (difference in the weight of the cuttings before and after saturation) or directly 
by NMR measurement. 

Figures 2a to 5a Figs. 2A to 5A show the result of experiments carried out from 
model cuttings, 1 to 2 mm in size, manufactured in the laboratory from rocks of known 
20 properties for which a curve Pc conventionally measured by centrifugation of a core is 
available. It can be seen that a good agreement is obtained with the reference curve at 
the level of the asymptotic part (low wetting fluid saturation). On the other hand, a big 
difference is observed for higher wetting fluid saturations. Results equivalent to those 
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obtained within the context of porosimetry measurements using mercury on cuttings are 
thus obtained, without pollution risks. 

The measured c apillary pressure curve Pc is however representative only on the 
asymptotic part. A reconstruction procedure is therefore necessary to evaluate the 
5 b e haviour behavior of the curve over the total saturation range. 

II) Measurement of the permeability k of the cuttings 

The method described in the aforementioned French patent application 02/02,242 is 
applied to measure the permeability of the cuttings. The cuttings are therefore dipped in 
a containment enclosure containing a viscous fluid. The enclosure is then 

10 communicat e d with coupled to a vessel containing the same fluid under pressure, so as 
to compress the gas trapped in the pores of the rock. According to a first embodiment, 
this communication period can be very short and followed, after a latency time, by the 
measurement of the pressure evolution in the enclosure. According to another 
embodiment, the communication period can be long enough to allow to obs e rv e 

15 observing and measure measurement of the variation of the volume actually absorbed 
by the cuttings. 

The evolution of the pressure or of the volume in the enclosure is then modelled 
from initial values selected for the physical parameters of the cuttings, and the values of 
the physical parameters of the cuttings are iteratively adjusted so that the modelled 
20 evolution best adjusts with the measured evolution of the physical parameter in the 
enclosure. 

This procedure gives excellent results. The permeability values k of the cuttings are 
totally in accordance with the reference measurements obtained from cores. 
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III) Reconstruction of the total capillary pressure curve 

In this third stage, w e synthe s ize the previous measurements are synthesized that 
have allowed to construct construction of the asymptotic part of the capillary pressure 
curve Pc and the measurements of permeability k using empirical relations r e put e d to 
5 model as well the physical parameters of rocks^veU. The following publications : 

- Thomeer, J.H.M. : "Introduction of a pore geometrical factor defined by the 
capillary pr e ssure curve Of A Pore Geometrical Factor Defined By The Capillary 
Pressure Curve ", Trans AIME, vol. March, pp.73-77, I960, and 

- Thomeer, J.H.M. : "Air p e rm e ability as a function of thre e pore network 
10 paramet e rs Permeabilitv As A Function Of Three Pore Network Parameters ", Trans 

AIME, vol. April, pp.809-814, 1983 

describe methods for evaluating the permeability from a capillary pressure curve. A 
Capillary capillary pressure curve Pc is modelled in the following form : 

P c =P d x exp( ) 

15 where : 

G i s a shape p arameter for taking account of accounting for the curvature o f the 
capillary pressure curve (related to the shape of the pore size distribution)^ 

- P d ]s_the displacement pressure extrapolated to S Hg equal to zero^; and 

V b oo is the percentage of volume occupied by the mercury at the end of the 
20 experiment at an infinite capillary pressure (equal to <|> x S Hg )^ 
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The three parameters of the model being related to the permeability by the 
following expression : 



In the following publication : 

5 - Swanson, B.F. : "A simpl e corr e lation between permeability and m e rcury capillary 
pres s ures Simple Correlation Between Permeability And Mercury Capillary Pressures ", 
JPT, vol. December, pp.2498-2504, 1981, 

the author proposes correlating the value of the permeability with the maximum value 
of ratio (Vb/P c ) on the mercury porosimetry curve. This particular point generally 

10 corresponds to the regime change that occurs at the end of the percolation regime, just 
before the significant capillary pressure increase. As regards the capacity of fluid flow 
in the porous medium, this point is particularly important because it represents the pore 
size for which the entire pore network is connected and which therefore controls the 
flow. The most general correlation provided by the author is given by the expression 

15 (Swanson 1981) : 



The following publication : 
- Kamath, J. : "Evaluation of accuracy of estimating air p e rmeability from m e rcury 



k =3.8068xG 



-1.334 



X 





2.005 




Of Accuracy Of Estimating Air Permeability From Mercury Injection 



20 Data ", SPE Formation evaluation, vol.7,4, pp.304-310, 1992 
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also relates to the comparative evaluation of the approaches using empirical correlation 
or physical models to determine the permeability value from a mercury porosimetry 
curve. The best agreement is obtained with a new correlation based on the 
aforementioned characteristic length by Swanson 1981, mentioned above : 



*=413xl££ if k<\mD 

k = 347x4*° ifk>\mD, 

L max being defined as follows : 



Anax - 



_ <*xlx(100-S r ) (pa = s w -s, 

A+i [r] 100-5, 



P e x(l + A) 



where : 

- X : exponent expressing the curvature of the capillary pressure curve (related to the 
10 shape of the pore size distribution), 

- P c : the displacement pressure extrapolated to S Hg equal to zero, and 
S r : residual saturation occupied by the wetting fluid (%). 

To parametrize the empirical capillary pressure curve Pc obtained from Thomeer's, 
Swanson's orKamath's approaches, it is constrained to adjust to the asymptotic part 

15 obtained by centrifugation during the first stage of the method. The entire curve is 
constrained using also the value of permeability k measured on cuttings during the 
second stage of the method, which is compared with the result of the empirical 
relations. The parameters of the capillary pressure Pc are then modified until both the 
measured asymptotic behaviour behavior and the permeability estimation are met, 

20 which allows to con s train constraining the capillary pressure curve over the total 
saturation range Sw. 
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The inversion process starts with a set of default parameters which allow te 
calibrat e calibration of the asymptotic behaviour behavior of the capillary pressure Pc 
with low water saturations. These parameters are then modified step by step (mainly the 
input pressure Pe or Pd and the shape factor X or G) so that the estimation of the 
5 permeability given by one of the previous relations is in good agreement with the 
permeability measurement obtained on cuttings while keeping a good agreement with 
the measurements of Pc at low water saturations. 

1.1 Results obtained 

Figures 2h. 3h. 1b and 5b 2B. 3B, 4B and 5B show the comparison between the 
10 curve Pc reconstructed according to the previous procedure and the reference curve Pc 
obtained on a core. It can be seen that, whatever the example considered, the 
reconstruction method allows to obtain obtaining a pertinent evolution of Pc over the 
total saturation range and in particular at high water saturations, whatever the 
permeability of the rock. 
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0.89 


0.38 


GDV1 


195 


0.28 


1.00 


0.41 


B7 


780 


0.25 
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0.26 



15 



0 



PATENT 



INSTITUT FRANQAIS DU PETROLE 

METHOD OF EVALUATING THE CAPILLARY PRESSURE CURVE 
OF ROCKS OF AN UNDERGROUND RESERVOIR 
FROM MEASUREMENTS ON DEBRIS 

Inventors : Roland LENORMAND, Patrick EGERMANN and Daniel LONGERON 

ABSTRACT 

Method A method of evaluating the capillary pressure curve of rocks of an 
underground reservoir from measurements on rock debris or fragments such as cuttings 
from the reservoir, over the total saturation range of these rocks, within a short period 
and at a low cost, from these measurements is disclosed . The method 

It essentially comprises measuring the permeability k of the debris, measuring the 
capillary pressure curve Pc as a function of the saturation of these fragments initially 
saturated with a fluid (brine for example) by subjecting them to centrifugation, and 
parametrizing a capillary pressure curve Pc satisfying empirical relations depending on 
adjustable parameters, constrained to adjust to an asymptotic part of the capillary curve 
measured by centrifugation, and to the value of permeability k measured on the cuttings, 
so as to obtain the whole of the capillary pressure curve. 

Applications .- include hydrocarbon reservoir evaluation for exampl e. 
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